Qutrit: entanglement dynamics in the finite qutrit chain in the consistent magnetic 

field 



O 



(N 



Oh 



c3 



X 



E. A. IvanchenkcQ 

Institute for Theoretical Physics, National Science Center "Institute of Physics and Technology" , 
1, Akademicheskaya str., 61108 Kharkov, Ukraine 
(Dated: June 14, 2011) 

Based on the Liouville-von Neumann equation, we obtain closed system of equations for the 
description of a qutrit or coupled qutrits in arbitrary time-dependent external magnetic field. The 
dependence of the dynamics on the initial states and magnetic field modulation is studied analytically 
and numerically. We compare the relative entanglement measure's dynamics in the bi-qutrit system 
with permutation particle symmetry. We find the magnetic field modulation which retains the 
entanglement in the system of two coupled qutrits. Analytical formulas for entanglement measures 
in the chain from 2 to 6 qutrits are presented. 



PACS numbers: 03.67.Bg Entanglement production and manipulation 
03.67.Mn Entanglement measures, witnesses, and other characterizations 
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I. INTRODUCTION 



Multi-level quantum systems are studied intensively, since they have wide applications. Some of the existent ana- 
lytical results for spin 1 [ll are derived in terms of the coherent vector The class of exact solutions for a three-level 
system is given in Rcf. [3[. The application of coupled multi- level systems in quantum devices is actively studied 
[j]. The study of these systems is topical in view of possible applications for useful work in microscopic systems [f|. 
Exact solutions for two uncoupled qutrits interacting with vacuum are obtained in Ref. @ . For the case of the qutrits 
O" 1 interacting with stochastic magnetic field exact solutions are obtained in Ref. Q • Exact solutions for coupled qutrits 
in magnetic field as far as we know were not found. 

The entanglement in multi-particle coupled systems is an important resource for many problems in quantum infor- 
■ mation science, but its quantitative value is difficult because of different types of entanglement. Multi-dimensional 
f — ' entangled states are interesting both for the study of the foundations of quantum mechanics and for the topicality 
0^ ■ of developing new protocols for quantum communication. For example, it was shown that for maximally entangled 
states of two quantum systems, the qudits break the local realism stronger than the qubits and that the entangled 
qudits are less influenced by the noise than the entangled qubits. Using entangled qutrits or qudits instead of qubits 
\& • is more protective from interception. From the practical point of view, it is clear that generating and saving the 
entanglement in the controlled manner is the primary problem for the realization of the quantum computers. The 
maximally entangled states are best suited for the protocols of quantum telcportation and quantum cryptography. 
The entanglement and the symmetry are the basic notions of the quantum mechanics. We study the dynamics of mul- 
tipartite systems, which are invariant at any subsystem permutation. The aim of our work is finding exact solutions 
for the dynamics of coupled qutrits interacting with alternating magnetic field as well as the comparative analysis of 
the entanglement measures in the chain of qutrits. 

The rest of the paper is organized as following. The Hamiltonian of the anisotropic qutrit in arbitrary alternating 
magnetic field is described in Sec. II. Then the system of equations for the description of the qutrit dynamics is 
derived in the Bloch vector representation. We introduce the consistent magnetic field, which describe entire class of 
field forms. In section III we derive the system of equations for the description of the dynamics of two coupled qutrits 
in the consistent field and find the analytical solution for the density matrix in the case of anisotropic interaction. 
Analytical formulas, which describe the entanglement in spin chains from 2 to 6 qutrits, are presented in Sec. IV. 
The results are demonstrated graphically in Sec. V for concrete parameters. The brief conclusions are given in Sec. 
VI. The auxiliary analytical results are presented in the Appendices. 
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II. QUTRIT 
A. Qutrit Hamiltonian 

We take the qutrit Hamiltonian (for the spin-1 particle) in the space of one qutrit C 3 in the basis |1 >= 

(1,0,0), |0 >= (0,1,0), | — 1 >= (0,0,1), in external magnetic field h = (hi, h 2 , h 3 ) with anisotropy, in the 
form 

H = fnSi + h 2 S 2 + h 3 S 3 + Q{Sl - 2/3E) + d{S 2 - S?), (1) 

where hi, h 2 , h 3 are the Cartesian components of the external magnetic field in the frequency units (we assume 
h = 1, Bohr magneton ps = 1); Si, S 2l S 3 are the spin-1 matrices (see Appendix A); E stands for the 3x3 unity 
matrix; Q, d are the anisotropy constants. When the constants Q, d are zeros, then the two Hamiltonian eigenvalues 
are symmetrically placed in respect to the zero level. 



B. Liouville-von Neumann equation 



The qutrit dynamics in the magnetic field we describe in the density matrix formalism with the Liouville-von 
Neumann equation 



idtp = [H, p], p(t = 0) = p . 



(2) 



It is convenient to rewrite Eq. ([2]) presenting the density matrix p in the decomposition with the full set of orthogonal 
Hermitian matrices C a @ (further the summation over Greek indices will be from to 8 and over the Latin ones 
from 1 to 8, see Appendix A) 
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(3) 



Since TrCi =0 for 1 < i < 8, then from the condition Tr p — Ro it follows that Rq = 1. And although the results are 
independent of the basis choice, in this basis the functions i?; = TrpC, have the concrete physical meaning [l(|. The 
values Ri,R 2 ,R 3 are the polarization vector Cartesian components; R4 is the two-quantum coherence contribution 
in R 2 ; i?5 is the one-quantum anti-phase coherence contribution in R 2 ; Rq is the contribution of the rotation between 
the phase and anti-phase one-quantum coherence; -R7 is the one-quantum anti-phase coherence contribution in Ri; 
i?8 is the two-quantum coherence contribution in R\ . 
Under the unitary evolution the length of the generalized Bloch vector 

(4) 

is conserved. The length of the generalized vector for pure states equals to Since id t p n = [H,p n ] (n = 

1,2,3,...), then under unitary evolution there is countable number of the conservation laws Tr p = Ci = 1, Tr p 2 = 
c 2 , . . . , from which only c 2 , c 3 are algebraically independent (Tlj . Additional quadric invariants of motion can be 
easily obtained after equating the matrix elements in defining the pure state. For example, two of these invariants, 
which follow from the expression (p 2 — p)i 3 = 0, have the form 



ti x - K 2 



Rl-R 2 - 2^|(1 - V2R 6 )R 8 = 0, R 5 R 7 - R X R 2 + ^=(4= - Re)Ri = 0. 



(5) 



For numerical calculations, these invariants control also the signs of the values Ri and thus the using of the invariants 
is useful when the analytical solutions are difficult to find. According to the Kelly-Hamilton theorem, the density 
matrix p satisfies to its characteristic equation 



P 3 -P 2 



2-b< 



p - detpE = 0. 



(G) 
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From equation (j6)) it follows that the density matrix determinant det p = (Trp 3 — Trp 2 )/3 + (2 — 6 2 )/18 is also the 
motion invariant. The Liouville-von Neumann equation in terms of the functions Ri takes the form of the closed 
system of 8 real differential first-order equations. This system of equations in the compact form can be written as 
following [ll|, Q3: 

dtRi = eijihiRj, (7) 

where eiji are the antisymmetrical structure constants, hj = 2(hi, fi3, 0, 0, -^=, 0, d) are the Hamiltonian components 
(P) in the basis C a (see Appendix A) . 

C. The consistent field 

Consider the qutrit dynamics in the alternating field of the form 

h(t) = (wicn(u;i|fc), u}±sii(ujt\k), ajodn(wt|fc)) , (8) 

where cn, sn, dn arc the Jacobi elliptic functions flU ] . Such field modulation under the changing of the elliptic 
modulus k from to 1 describes the whole class of field forms from trigonometric (cn(wt|0) = cosut, sn(wi|0) = 
sinwi, dn(wi|0) = 1 ) [13] to the exponentially impulse ones (cn(w£|l) = -^-r, sn(cji|l) = thwt, dn(o;£|l) = ^jj) 
[l5j . The elliptic functions cn(wi|fc) and sn(ujt\k) have the real period ^jf. while the function dn(u>t\k) has the two 
times smaller period. Here K is the full elliptic integral of the first kind 0]. In other words, even though the field is 
periodic with common real period , but as we can see, the frequency of the longitudinal field amplitude modulation 
is two times higher than the one of the transverse field. Such field we call consistent. 

Let us make use of the substitution p = a^ 1 ra\. Then we obtain the equation for the matrix r in the form 

idtr = [aiHa^ 1 - iaidtia^ 1 ), r] (9) 

with the matrix 




a x = | 

where f(ujt\k) = cn(ojt|fc) + isn(a;t|fc). Since 

aiSic*! 1 = Sicn(uit\k) — 52sn(wi|fc), ai^aj -1 = Sisn(u;i|fc) + S2cn(uit\k) , aiS^a^ 1 = S3, (11) 
then the equation for the matrix r without taking into account the anisotropy can be written as following 

idtr = [toiSi + S dn(wi|fc)S l 3, r], r(t = 0) = po, S = ojq — uj. (12) 
At k = equation (|12p describes the dynamics of the qutrit in the circularly polarized field [3, [l6|, [13] • The exact 



solutions of this equation are known and at some initial conditions the explicit formulas are given in Ref. [18J . At 
exact resonance, w = ujq it is straightforward to present Eq. ^ in the deformed field (k 7^ 0) ([5]) for the given initial 
condition p = p : 

p(t) ^a^e-^^poe^^aL (13) 

Explicit solutions for some specific initial conditions are given in the Appendix B, Eqs. (|3"9"|) - (HH) . From the 
explicit exact solutions in the deformed field at resonance <5 = one can see that the populations and the transition 
probabilities do not depend on the field deformation (it is independent of the k modulus). 
Consider the solution of Eq. (|12[) far from the resonance in the form of 6 power expansion 

r(t) =A°\t)+r^(t) + --- . (14) 

Then we put the expansion (fT4|) in Eq. (fT2j) and equate the same degree terms. As the result we obtain the system 
of equations for finding r^'(£): 



ifitrW =oj 1 [S 1 ,A% (15a) 
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idtr® =u) 1 [Si,rW}+6dn(u;t\k)[S 3 ,r( l -% 1 = 1,2, .... (15b) 

We multiply Eq. (|15[) to the left by the matrix e luJ1 1 and to the right by the matrix e ~ 4a;itSl for formation of the 
integrating multiplier (l9j . Now finding the terms in the scries (|14[) is defined by the previous ones r^ _1 ) as 
following 

r w (t) = -i$ f dt'e iul(t '- t)Sl dn(ujt'\k)[S3,A l - 1 \t')}e- iu ^ t '- t)Sl . (16) 
Jo 



III. BI-QUTRIT 

In the space C 3 <g> C 3 the two-qutrit density matrix can be written in the Bloch representation 

Q = -^RapCa <g> Cfj, Roo = 1, g(t = 0) = g , (17) 

where ® denotes the direct product. The functions R m o,Rom characterise the individual qutrits and functions R mn 

characterise their correlations. The length of the generalized Bloch vector \J R^p ~ ^ ^ or P ure states equals 2\/2. 

Consider the Hamiltonian of the system of two qutrits with anisotropic and exchange interaction in magnetic field in 
the following form 

H 2 = (t~3 + Q(S% - 2/3E) + d(Sf - Sf )) ®£ + 
E®(h~S" + Q{Sl-2/?>E) + d{Sl- Sl)) + JSi®Si = ^h a[i C a ®C fjl (18) 

where h and h are the magnetic field vectors in frequency units, which operate on the first and the second qudits 
respectively, and J is the constant of isotropic exchange interaction. 

The system of equations for two qutrits takes the real form in terms of the functions R m o, Rom, Rmn as the closed 
system of 80 differential equations [l2j , supplemented by the initial conditions 

[2 [2 
dtRmO = Y ^epimihpoRio + hpiRu), d t Ro m = y -^ e pim{hopRoi + hipRn), (19a) 



dfRr, 



-*pim 



— (h pn RiQ + h p0 Ri. 



QrlnhprRi, 



+ e 



pi II 



(19b) 



where by definition 



Tr P C a ®C[} = -R at 



(20) 



and h 



V6(t,0,0,%,0,d), hop 



p0 



s/3> 



s/6(h,0,0,^=,0,d), h n 



'33 



2 J are the Hamiltonian expansion 



coefficients in the basis C a ® Cp (other coefficients equal to zero). In equations (|19[) Latin indices m, n take the values 
from 1 to 8. Numerical values for the structure constants e P i m , g r im are given in Appendix A. 
The energy of the coupled qutrits in terms of the correlation functions has the following form 



E(t) 



3 



(h p oRpQ + hopRop + huRitj 



(21) 



i=l 



We study the dynamics of two qutrits in the magnetic field h = (u>icn(u>t\k)), uisn(uit\k), o;odn(wt|fc), h 
(z&icn(mt\k) , vjisn(ujt\k)i vuo<hi(ujt\k)) at the anisotropy constants equal to 0. We transform the matrix density 
q = Q' 2 " 1 r2a2 with the matrix ui = a.\® ot\. Then equation for the matrix r-i takes the form idtr^ = [H, with the 
transformed Hamiltonian 
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Since £> = chi(iL)t\k)\k = o = 1, then the transformed Hamiltonian H does not depend on time and the solution for the 
density matrix in the circularly polarized field has the form 



g(t) = a 2 



-iHt iHt I 

£>oe a 2 \k=o- 



(22) 



In the consistent held at resonance lu = zdq = uiq = h at equal vj\ = uj\ the Hamiltonian eigenvalues equal to 
— 2J, — J, J, J — 2oji, — J — wi, J — u>i, — J + uj\, J + wi, J + 2lj\. This allows to hnd the exact solution in the closed 
form for any initial condition, since the matrix exponent e tHt in this case can be calculated analytically. 
For larger number of the qutrits with pairwise isotropic interaction, the generalization is evident. In the case of 
interaction of qudits with different dimensionality, the reduction of the original system to the system with constant 
coefficients can be done by choosing, for example, the transformation matrix for spin-3/2 and spin-2 in the form 



dia g (/ 3 / 2 , f\ r lj \ r 3/2 )^dia g (/ 2 , /, i, r\ r 2 )- 



(23) 



However, the Hamiltonian eigenvalues cannot be found in the simple analytic form because of the lowering the system 
symmetry. 



IV. ENTANGLEMENT IN THE QUTRITS 

A. Entanglement in the bi-qutrit 

For the initial maximally entangled state, which is symmetrical at the particle permutation, 

I 1 

\i> >= -7= E N > ®l* >. ( 24 ) 

v i— — 1 

in the consistent field at the resonance w = wq = uiq = h at equal w\ = wi, the exact solution for the correlation 
functions is given in Appendix C. The correlation functions have the property R a p = Rp a , i.e. the symmetry is 
conserved during the evolution, since the initial state and Hamiltonian are symmetric in respect to the particle 
permutation. 

Given the exact solution, one can find the negative eigenvalues of the partly transposed matrix g pt = (T <E> E)g (here 
T denotes the transposition): 

ei = e 2 = -^\/69 + 28cos3Jt - 16cos6Jt, e 3 = -— (5 + 4cos3J/j) . (25) 

The absolute value of the sum of these eigenvalues myw — |ei + £i + £3 1 defines the entanglement measure (negativity) 
between the qutrits [2(| ■ 

The entanglement between the qutrits can be described quantitatively with the measure [21 1 



msM = \j\{Rij - R l0 Roj) 2 - (26) 

This measure equals to for the separable state and to 1 for the maximally entangled state, and it is applicable for 
both pure and mixed states. 

That is why for the maximally entangled initial state of two qutrits, the entanglement in the consistent field is defined 
by the formulae with the found solution for the density matrix 

m SM = 1 V4457 + 2776 cos 3 Jt - 632 cos 6 Jt - 56 cos 9 Jt + 16 cos 12 Jt. (27) 
V6561 
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This measure is numerically equivalent to the measure myw 12QL l22[, which is defined by the absolute value of the 
sum of the negative eigenvalues (|25p of the partly transposed matrix. 

According to the definition [23| for 2-qutrit pure state, the entanglement measure equals to 



2 

1 , 

m 



|l>, ( 28 ) 



2 

i=i 



where Si = — Tr pi log 3 pi is the reduced von Neumann entropy, the index i numerates the particles, i.e. the other 
particle are traced out. 

Since the qutrit reduced matrix eigenvalues equal to A x = A 2 = ^(5 + 4cos3Jt), A 3 = ^(17 — 8cos3Jt), then the 
entanglement measure in the bi-qutrit takes the form 



^ A, log 3 A,. (29) 



Normalized by the unity, the measure I-concurrence, which is easy to calculate, is defined by the formulae [24 



mi = ^— J2(l -Trp 2 ) = -V57 + 32 cos3 Jt - 8 cos 6 Jt, (30) 
2 * 9 



where pi = -^=C a R a o is the reduced qutrit matrix. 

The timc-dcpcndcncc of the measure msM for the symmetrical initial state 



12 . 

in 



takes the form 



SM 



1 



V209952 



V102679 + 19136 cos 3 Jt + 29312 cos 6 Jt - 1024 cos 9 Jt + 800 cos 12 Jt; (32) 



at t = this measure equals to ^23/32. 

The measures myw, msM, i]2, nil, mg M do not depend on the parameters of the consistent field, sign of the exchange 
constant at zero anisotropy parameters. It should be noted that the Woofers entanglement measure (concurrence) in 
the system of two qubits with the isotropic interaction in the circularly polarized field at resonance is also independent 
of the alternating field amplitude [25j . but depends on the exchange constant magnitude and the initial conditions 
only. 

At zero external field the entanglement measure (|24[) takes the analytic form at equal non-zero anisotropy parameters 
Q=d=d=Q 



msM{Q) 



(9J 2 + 8QJ+16Q 2 ) 2 \ 



^2 q k cos (ky/9 J 2 + 8Q J + 16Q 2 1) , (33) 



k=0 



where q = 4457 J 8 + 11616Q J 7 + 47392Q 2 J 6 + 85888Q 3 J 5 + 163072Q 4 J 4 + 194560Q 5 J 3 + 221184Q 6 J 2 + 131072Q 7 J + 
65536Q 8 ; q x = 8 J 2 (J + 2Q) 2 (347 J 4 + 518QJ 3 + 1440Q 2 J 2 + 1504Q 3 J + 1024Q 4 ) ; 
q 2 = -8J 2 (J + 2Q) 2 (79 J 4 + 76QJ 3 + 320Q 2 J 2 + 448g 3 J + 256Q 4 ); 
q 3 = -8J 3 (7J - 4Q)( J + 2Q) 3 ( J + 4Q), q 4 = 16J 4 (J + 2Q) 4 . 

B. Entanglement in the chain of qutrits 

We consider the Hamiltonian of the chain of N qutrits with the pairwise isotropic interaction in the magnetic field 
at in the following form 

N-l N-2 

H N = ^(^ ® E ® • • • ® E+J~t ® ~t ®E <g> • • • <g> E), (34) 
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FIG. 1: The time-averaged populations for the initial pure state | — 1 > versus the normalized Larmor frequency ujq/lj at the 
parameters k = 0.85 (solid line), k — 0.2 (dashed line), d — Q — 0, Ui = 1/3, oj = 1 (/ shows the upper level |1 > population; 
77 shows the middle level |0 > population). 



where the summation is over different possible positions of ~is in the direct products. Because the maximally entangled 
state of N qutrits 

i^>=7fEN>^ ( 35 ) 

and the Hamiltonian (|34p have the permutation symmetry, it follows that the density matrix of N qutrits has the 
symmetric correlation functions. The length of the generalized Bloch vector for pure states equals s/3 N — 1. 
The entanglement measures for the many-particle multi-level quantum systems are not studied enough and difficult 
to calculate in the analytic form, that is why we will present only analytic formulas for the entropy measure t]n p3j . 
which is defined by the eigenvalues of the reduced one-particle matrices for each qutrit. As the result of the mentioned 
symmetry, the reduced matrices are equal to each other. Therefore the entanglement measure for N qutrits is defined 
by the formulae 

3 

Vn = n log 3 n. (36) 

i=l 

The eigenvalues of the reduced matrices for 3, 4, 5, and 6 qutrits are presented in the table below 

N \ n n = r 2 r 3 

g 29-4 cos 5 Jt 17+8 cos 5Jt 

a 905—98 cos 3 Jt — 72 cos 7 Jt 395+196 cos 3 7 Pt+144 cos 7 Jt /Q7\ 

2205 2205 V 5 ' > 

r 16919-1944 cos 5 Jf -800 cos 9 Jt 8687+3888 cos 5 Jt+1600 cos 9 Jt 



21977-1694 cos 3 Jt-'W^cos 7Jt-560 cos 11 Jt 9407+3388 cos 3Jt+3l^ 2 cos 7Jt+1120 cos 11 Jt 
53361 53361 



The measures 773, 774, 775, tjq do not depend on sign of the exchange constant like the measure 772- 



V. NUMERICAL RESULTS 



In Fig. [T] we present the populations of the upper and middle levels in the qutrit averaged over the time interval 
t — y 00: P+ = i J r dt + -^i?3(i) + 575 -ReO)) ,P° = 7 Jo dt{\ - \V2R & {t)) in dependence on the normalized 
Larmor frequency luq/uj. The population of the upper level in qutrit coincides in form with the upper level occupation 
in a two-level system [l9j |. i.e. this demonstrates the magnetic resonance position stabilization and the presence of 
the parametric resonances. 

In Fig. [2] we note the considerable suppression of the qutrit spin oscillations S y = cn(wi|fc) sinwiTj and S z — — cosuiit 
by the environment (fluctuator) in the case of the resonance u) = uq, w = 0. 
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120 t 240 



FIG. 2: Dynamics of the spin vector components S y , S z for the initial pure state | — 1 > (dashed lines) in the circularly 
polarized field with the parameters: k — 0, u>i = 0.02, u> = cjq = 1, d — Q = 0. Solid lines demonstrate the deformation of the 
spin components due to the influence of the second spin (the fluctuator) with J = 0.1 for the initial pure state | — 1 > ®| — 1 >. 




FIG. 3: Disentanglement dynamics of the initially maximally entangled state in the bi-qutrit: in the zero external field with 
equal anisotropy constants Q = d = d = Q = 0.02507, J = —0.1 (curve 1) and for J = 0.1 (curve 2); in the consistent field 
the curve 3 (thick line) demonstrates complete coincidence of the measures mvw and msM at J=0.1 and zero anisotropy 
constants; the curve 4 demonstrates the entropy measure 772; /-concurrence is presented by the curve 5 at J— 0.1. 

The bi-qutrit energy ([21]) in the consistent field at isotropic interaction in the case of the solution (|45|) is constant 
and equal to | J. 

Although the analytic expressions for the measures in a bi-qutrit mywi tusm are different, but the numerical values 
are practically identical. Maximal deviation in the rectangle (1 > J > 0.01) x (100 > t > 0) equals 0.014, where x 
denotes the Cartesian product. 

Measures 772 and mj qualitatively coincide with the measures mvw, mSM- 

We have found that the anysotropy of the qutrits disentangles them, namely the entanglement is decreased down to 
0.0010 (see graphs 1 and 2 in Fig. ©. 

In the constant longitudinal field at = — w = (0, 0, cjq) (the bi-qutrit Hamiltonian eigenvalues are equal to 
J, J, Xi,X2, X3, —p, —p,p,p, where Xi,X2,X3 are the roots of the equation x 3 + 2x 2 J — p 2 x — 2J 3 = 0, p = \J J 2 + uj 2 ) 
the Hamiltonian contains the antisymmetric part, thus it follows that the density matrix for the initial symmetric state 
will not be symmetric because of the breaking the symmetry of the particle permutations. The analytic solution is cum- 
bersome. In the constant longitudinal impulse field "d* = = (0, 0, 2 (6((t- 17) (t- 60)) + (9((40-i)(57-i)(t-60)))) 
the entanglement dynamics is blocked at luq J (Fig. |4]). This points to the possibility to control the entanglement. 
In Fig. [5] we present the comparative dynamics of the entropy entanglement measure for 2 to 6 qutrits. The disen- 
tanglement dynamics of the measures 773 , 774 , 775 , rje is similar to the one in the case of two qutrits, but with smaller 
oscillation amplitude, i.e. larger number of the qutrits disentangles less than two qutrits (0.889 < 773 < 1). 

VI. CONCLUSION 

We have shown that the time-averaged upper-level occupation probability for the qutrit in the consistent field in 
dependence on the normalized Larmor frequency loq/uj coincides in form with the upper-level occupation probability 
in the two- level system and the parametric resonances appear (Fig. [1]). In the qutrit coupled to another qutrit 
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FIG. 4: Disentanglement of the maximally entangled state (|24[) (solid line) in the impulse field idi 
2, J = 0.f78. The dashed curve presents the evolution in zero external field. 



0, Cl>0 



-zuo 




FIG. 5: Disentanglement of the maximally entangled state ()28[) in the chain of 2,3,4,5 and 6 qutrits with J = 0.1. 



(fluctuator), the spin oscillations are essentially suppressed. 

The comparative analysis of the bi-qutrit entanglement measures on the base of the analytic solution for the density 
matrix demonstrates that, in spite of the different approaches to the derivation of the formulas for the entanglement, 
all the formulas give quite close results (Fig. [3]), and the measures myw and tusm are practically equal. This is in 
accordance with the general results for the entanglement in the systems with the permutational symmetry [22[ ■ 
The analytical formulas for the entanglement measures 773 , 774 , 775 , 776 are similar to the disentanglement measure for 
two qutrits 772 , but with numerically smaller oscillation amplitude, i.e. the larger number of the qutrits disentangles 
less than two qutrits. 
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VII. APPENDIX A 



The matrix representation of the full set of Hermitian orthogonal operators for spin-1 has the form 






C% — S2 — — p= 
V2 



/CI -1 
1 
1 

V 








(38a) 



10 



c 4 



0-1 



1 



7 C5 




C 6 = V3(Sf - 2/3E) 



V3 



/ 1 

0-2 
1 



(38b) 




C_ q2 q2 
8 — 1 — 2 — 




C = \/-E, 



(38c) 



where E is the unity 3x3 matrix. These matrices have the property of the trace equal to zero Tr C a = and 
orthogonality Tr C a Cb = 25 a b, 1 < a, b < 8. The connection between the basis C a and the Gcll-Mann basis A a is the 
following: 

C\ = ~y|(^i + ^6)7 C2 = (A2 + A 7 ), C 3 = ^A 3 + -2^A 8 ,C 4 = A 5 , C 5 = -^5(^2 — A 7 ), C 6 = ^A 3 — ^Ag, = 
— Ag), Cg = A4. 

Non-zero structure constants e a bc (dabc) antisymmetric (symmetric) in respect to the permutation of any pair of 



indices for the commutators [C a ,Cb] = 2ie a bcC c (anticommutators {C a ,Cb} = ^E8 a b + 2g a b c C c ) are respectively 

ii 



equal according to the dchnitions e a bc = ji Tr [C a , Cb]C c : 
e 



123 — e 147 — ei58 — — e-245 — ^278 — ^357 — 5, e 15 g — e 2 67 — — j ^348 — — lj 9abc — 4 Tr {C a , Cb}C c : 5336 — .9446 
1 1 
—3666 — 3688 — TyJ: 3116 — 3226 — 3556 — 3677 — — 3118 — 3124 — 3137 — — 3228 — 3235 — ~ 3457 — 3558 — — 3778 

i. Hence, the product of the generators is equal to C a Cb = ^ES a b + (g a bc + *e a b c )C c . 



VIII. APPENDIX B 

For the initial state | — 1 > at non-zero detuning 5 = ujq — lu the density matrix elements in the circularly polarized 
field have the form 



fit 



\/2w? , fit 



fit , . fit\ 

%UJt 1 ifi cos \-6 sin — I 



2 J 



(39) 



Pl3 



2 ®t „-2iu>t 



— — — sm" — e 

2fi 4 2 



2iSn sin fit + (25 2 +w 2 ) cos fit) , 



(40) 



w?sin 2 ^ . „ „ ,. 
P22 = ^ 4 2 (2<5 2 + w 2 (l + cos fit)) , p 23 



V2fi 



(2<5 2 + w 2 (l + cos fit)) ( (5 sin 2 y +i- sin fit 



1 2 
P33 = "77^4 (2<5 2 +^ 2 (1 + COS fit)) , Plk = P * u , 



(42) 



where fi = \/ u\ + S 2 is the Rabi frequency. 

For the initial doubly stochastic state -j^-fl - 1 > +|0 > +|1 >) and for the states |0 >, || — 1 > +-^|0 > +||1 > 
at exact resonance 6 = in the consistent field, the density matrices are respectively equal 



^ (cos 2wii + 3) i/" 1 (iV2sin2a;it + 4) 

^/(4-« v / 2sin2a;it) f (3 - cos 2wit) 
^/ 2 (cos2a;it + 3) ±f (iy/2 sin 2uiit + 4) 



±f~ 2 (cos2wi* + 3) 
j^f- 1 (4- i-v/2sin2wit) 
fk (cos2wit + 3) 



12 



(43) 



/ I 



sin^t - ^" lBi ° 2 ^ t l/- 2 sin 2 ^t\ / ^(5-cos2^t) i/^sin 2 ^ 



i/ sin 2ui£ 
2\/2 

y i/ 2 sin 2 ^ 



2%/2 
COS 2 Wit 

sin 2ui£ 

2~72 



2^2 

i sin 2 wit 



/ 



if sin 2u\t 
8\75 

^ |/ 2 sin 2 wit 



8%/2 8 _ 

|(cos2wit + 3) ir \% 2uilt 

if sin 2a^ij 

8^2 



j£ (5 — cos2wit) , 



(44) 



For both the initial middle level and the doubly stochastic pure initial state, the populations of the upper and bottom 
levels are equal. [HI . This property is fulfilled for the mixed state as well. 
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IX. APPENDIX C 



The exact solution for the correlation functions of the initial state (|24p . which is symmetric under the particle 
permutation, in the consistent field at resonance w = Wq = ujq = h and equal W\ = lo\ takes the form 

8 [2 2 2 3Jt 4/2 2 3Jt 

Ro,i = Ro,2 = Rq,3 = 0, R0.4 = - y g cos wrf cnu snu sin i£ ,5 = — - y - cmt sin -^-sm2wit, 

2 /t: ,„ „ x . 2 3 „ 4/2 . „ . 9 3 Jt 

R o e = — v 2 (3 cos 2wit — 1) sin , /t 0i 7 = — y - snu sin2wit sm , 

9 2 3 V 3 2 

4 !~2 3 Jt 

R0.8 = g y g cos 2 wii (1 - 2sn 2 u) sin 2 — , (45a) 

fli! = — (16 + 12(cos 3 Jt + 2) (cn 2 u- sn 2 u) cos 2 wit + 2 cos 3 - 3 cos (3 J - 2wi) t 
36 v v 7 

2 

— 12 cos2wit — 3 cos(3 J + 2uii)t) , R\^ = - (cos 3 Jt + 2)cnusnu cos 2 wit, 

3 

Rx,s — — (cos3Jt + 2)snu sin2c<Ji£, -R1.4 = —snu sin3Ji sin2wi£, 
3 3 

Ri 5 = - (2 (cn 2 u — sn 2 u) cos 2 wit + 3cos2wit — l) sin 3 Jt, R\ % = —= cnu sin 3 Jt sin 2wit, 
6 V3 

2 . 1 

i?i t 7 = — — cos 2 wit cnu snu sin 3 Jt, i?i,8 = —7=R\fi, (45b) 

3 V3 

i?22 = — (6(cos 3 Jt + 2) (sn 2 u — cn 2 u) cos 2 w^t + cos 3 Jt — 3(cos 3 Jt + 2) cos 2w^ + 8) , 
18 

i?23 = --(cos3Jt + 2)cnu sin2wit, i? 2 4 = -4=i?i6, = —R17, R26 — V%Ru, 
3 V3 

i?27 = - (2 (cn 2 u — sn 2 u) cos 2 w x i — 3 cos 2wit + l) sin 3 Jt, i?28 = i=R26, (45c) 

6 V3 

i?33 = — (-2cos3Jt + 3cos(3J-2wi)i+12cos2wit + 3cos(3J + 2wi)t + 2), 
18 
2 

-R34 = — - cos 2 wit (cn 2 u — sn 2 u) sin 3 Jt, i? 35 = — -R14, -R36 = 0, 

1 4 

i?37 = — - cnu sin 3 Jt sin 2wit, i?38 = - cos 2 w^ cnu snu sin 3 Jt, (45d) 

i? 44 = — (-72 (l - 8cn 2 usn 2 u) cos 4 u)\t + 8 cos 3 Jt - 12(2 cos 3 Jt + 1) cos 2wif + 9 cos 4wif + 19) , 

i?45 = — snu (24 (sn 2 u — 3cn 2 u) sinwit cos 3 wit + 2(2cos3Jt + 1) sin2wit — 3sin4wit) , 
2 

i?46 = p cos 2 w x t (2 cos 3 Jt — 9 cos 2wi t + 7) cnu snu, 

3v3 

i?47 = — cnu (—24 (cn 2 u — 3sn 2 u) sin wit cos 3 Wit — 2(2 cos 3 Jt + 1) sin 2wit + 3 sin4wit) , 

i?48 = 4 cos wit cnu snu (cn 2 u — sn 2 u) , (45e) 

i?55 = — (6 (cos 3 Jt + 6 cos 2wit — 4) (sn 2 u — cn 2 u) cos 2 wit — cos 3 Jt + 3(cos 3 Jt + 2) cos 2wif 
18 

— 9cos4wit + 1) , i?56 = — -j= cnu ( 4sin 2 sin2w!t — 9sin4wit J , 

6V3 V 2 / 

i?57 = — (2 cos 3 Jt + cos (3 J — 2wi) t + 4cos2wit + 6cos4wit + cos (3 J + 2wi) t — 2) cnu snu, 

i?58 = — cnu (—24 (cn 2 u — 3sn 2 u) sinwit cos 3 wit + 2(2 cos 3 Jt + 1) sin2wit — 3 sin4wif) , (45f) 
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i? 66 = — (-4 cos 3 Jt + 6 cos (3 J - 2uj x ) t - 12 cos 2u\t + 27 cos 4u>if + 6 cos (3 J + 2wi) i + 13) , 
36 

i? 6 7 = — p sn-u (2(cos3Jt — 1) sin2wit + 9sin4wi<) , 
6y3 

i?68 = — 7= cos 2 wii (— 2 cos 3 + 9 cos 2wi t — 7) (cn 2 w — sn 2 u) , (45g) 
3v 3 

i? 77 = ^ ((cos3Jt - 1) (3cn 2 u - 3sii 2 m - l) + 9cos4o;it (cn 2 u - sn 2 w - l) 

+3(cos3J£ + 2) cos2wii (cn 2 u — sn 2 it + l)) , 
i?78 = — sou (6 (3cn 2 u — sn 2 u) cos 2 uj\t + 2 cos 3 Jt — 3 cos 2u)it + l) sin2u;ii, (45h) 

i?88 = (72 (l - 8cn 2 usn 2 w) cos 4 u;ii + 8cos3Jt - 12(2 cos 3 Jt + 1) cos 2oj x t + 9 cos4wit + 19) , (45i) 
where u = (ht\k). It is straightforward to find the analytic solution for larger number of qutrits at the same conditions. 
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